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ABSTRACT 

The ratio of the period of the fundamental mode to that of the first overtone of kink oscillations, from here on the ‘period ratio’, is 
a seismology tool that can be used to infer information about the spatial variation of density along solar magnetic flux tubes. The 
period ratio is 2 in longitudinally homogeneous thin tubes, but it differs from 2 due to longitudinal inhomogeneity. In this paper 
we investigate the period ratio in longitudinally inhomogeneous prominence threads and explore its implications for prominence 
seismology. We numerically solve the two-dimensional eigenvalue problem of kink oscillations in a model of a prominence thread. 
We take into account three nonuniform density profiles along the thread. In agreement with previous works that used simple piecewise 
constant density profiles, we find that the period ratio is larger than 2 in prominence threads. When the ratio of the central density 
to that at the footpoints is fixed, the period ratio depends strongly on the form of the density profile along the thread. The more 
concentrated the dense prominence plasma near the center of the tube, the larger the period ratio. However, the period ratio is found 
to be independent of the specific density profile when the spatially averaged density in the thread is the same for all the profiles. An 
empirical fit of the dependence of the period ratio on the average density is given and its use for prominence seismology is discussed. 
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1. Introduction 

High resolution observations reveal that solar prominences are 
formed by a myriad of long and thin sub-structures usually called 
threads or fibrils. In observations of prominences above the limb, 
vertical threads are commonly seen in quiescent prominences 
(e.g., iBereer et al.l 120081) . while horizontal thre ads are usually 


observed in active region prominences (e.g.,JOkamotoetjiI 


2007). However, it has been suggested (e.g., Schmieder et al. 


20101) that vertical threads might actually be a pile up of horizon¬ 
tal threads which appear as vertical structures when projected on 
the plane of the sky. Along this line of thought, threads supported 
in quasi-horizontal magnetic fields seem to be more consis¬ 
tent with observations of filament threads on the solar disk (see 
Linll201 lib with determ inations of the magnetic field orientation 
in pro minences (e.g., Casin i et aD 120031: lOrozco Suarez et akl 
2014a), and with equilibrium mode ls of prominences b ased on 
slightly dipped magnetic fields (e.g., Terrad as et al.l2013l) . In this 
work we use the term ‘thread’ to refer to prominence horizontal 
fine structures. 

Transverse oscillations of prominence threads with periods 
roughly b etween 1 and 20 minutes have been frequently reported 
(see, e .gJOkamoto et alJ2007t Lin et al.lf2007l 120091: iNing et akl 


2009; lOrozco Suarez et aL 2014bl) . From the theoretical point 
of view, the oscillations are usually interpreted as magnetohy 
drodynamic (MHD) kink modes of the flux tube tha t supports 
the prominence thread (see. e.g..lTbr radas etaL 2008^. On_etaI 


2009; 


2011 


SoleretaD 2010t lArregui et al. 2011 ; Soler & Goossens 


Soler et aHl2012 ). Both observational and theoretical 


pects of prominence thread oscillations have been reviewed by 
lArregui et al.l d2012bi) . 

Prominence seismology relies on the comparison of ob¬ 
served with predicted properties of prominence oscillations 
(see iBallesten l2014i) . The predicted properties are based on 
theoretical models. The aim is to indirectly infer information 
about th e plasma and/or the magnetic field in prominences 
(see, e.g., Terradas et al. 1120081 iLin et al. 1120091 ISoler et akl feo 10; 


as- 


lArregui et al.ll2012a ). In this direction, the ratio of the period 
of the fundamental longitudinal mode to period of the first 
longitudinal overtone of thread kink oscillations, from here 
on the ‘period ratio’, can be used to obtain information about 
the spatial variation of density along the threads. Standing 
MHD waves on magnetic flux tubes can be characterised by the 
number of nodes in their eigenfunctions, and this classification 
is independent of the velocity polarisation of the wave. Modes 
can have different number of nodes in the longitudinal (axial) 
part of the eigenfunctions as well as in the radial part of the 
eigenfunctions. The terms fundamental longitudinal mode and 
first longitudinal overtone refer to waves that have no nodes and 
have only one node, respectively, in the longitudinal (axial) part 
of the eigenfunction. In the present work we study the funda¬ 
mental longitudinal mode and the first longitudinal overtone of 
transverse kink waves. The period ratio is 2 in longitudinally 
homogeneous thin tubes, but i t differs from 2 due to longitu¬ 
dinal inhomogeneity dAndries et al.l l2005bllah . In the context 
of coronal loop transverse oscillations, the use of the period 
ratio as a seismology tool was first proposed by lAndries et al.l 
(2005a ) and has been exploited in a number of subsequent 
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2006; 

Donnellv et al. 

2007; Van Doorsselaere et al. [2007^ 

Verth & Erdelvi 2008 

; Ruderman et al. 2008; Arregui et al. 


among others). An extensive review of the use of the 
period ratio for coronal loop seismology can be found in 
lAndries et al.l (I2009I) . 

Although no reliable simultaneous observations of the two 
periods in prominence threads are curr ently available, numerica l 
simulations of thread oscillations by ISoler & Goossensl ( 201 1 ) 
indicate that the fundamental mode and the first overtone are 
the modes that are more easily excited in threads due to ex¬ 
ternal disturbances. It is reasonable to expect that the two pe¬ 
riods will be eventually detected in prominence threads by fu¬ 
ture high-resolution observations. Due to its great potentiality 
for seismology, a firm theoretical basis on the value of the pe¬ 
riod ratio in prominence threads is needed. The use of the period 
ratio for prominence seismology was first studied bv lDfaz et al.l 
( 20id ). They considered a prominence thread model composed 
of a thin magnetic flux tube filled with a dense blob of plasma 
that occupied a small part of the tube only, while the rest of the 
tube was evacuated, i.e., it was occupied by a much less dense 
plasma. The density in both the dense blob and in the evacu¬ 
ated part was homogeneous, and an abrupt jump between the 
two densities was assumed. Similar piecewise constant models 
have been used in previous works as, e.g., j pf 

ikld2C 


, e.g., iDIazet al.l 
Dvmova & Ruderman ( 20051): Terradas et al. (120081): Soler et alJ 


d200l: 


d2010l) : Soler & Goossens ( 2011 ): Soler et al. ( 20121) . while a 
continuous but thin transition between the dense blob and the 
e vacu ated part of the tube was considered by Arregui et al 
( 2011 ). The simple piecewise constant model used by Diaz et al 


(2Q10) allowed these authors to derive analytic expressions for 
the periods of the fundamental mode and the first overtone. They 
found that, contrary to the case of coronal loops where the period 
ratio is smaller than 2 due to stratification, in prominence threads 
the period ratio is larger than 2. The reason for this opposite re¬ 
sult is that in prominence threads the densest plasma is located 
near the tube centre, while in stratified coronal loops the pl asma 
is den ser near the footpoints. Subsequently, ISolei^&GoqssgnsI 
d201 ll) added longitudinal flow to the model of Diaz et alJ (12010) 
and f ound corrections to their results (see also Erdelvi et alJ 
[20f4l) . 


The purpose of this paper is to advance the study of the 
period ratio of kink oscillations in inhomogeneous prominence 
threads and to explore its implications for prominence seismol¬ 
ogy. Here we go beyond the piecewise constant density profile 
used by iDlaz et aH (1201 Oh and consider a more realistic contin¬ 
uous variation of density along the thread. Assuming that the 
thread formation is due to condensati on of cool plasma via ther¬ 
mal instability (see, e.g.. iLuna et al.l [2012b . the density profile 
along the threads is likely to depend on the energy balance be¬ 
tween footpoint heating, radiative losses, and thermal conduc¬ 
tion. This balance naturally leads to a spatially inhomogeneous 
density along the threads. Here, we evaluate the impact of the 
form of the longitudinal density profile on the period ratio. To 
do so, we numerically solve the full two-dimensional (2D) eigen¬ 
value problem of kink oscillations in a longitudinally inhomoge¬ 
neous prominence thread model. 

This paper is organized as follows. Section [2] contains the 
description of the prominence thread model and the numerical 
method used to solve the eigenvalue problem of kink oscilla¬ 
tions. The results of the computations are given in Section [3 
The implications of the findings of this paper for prominence 
seismology are discussed in Section [4] Finally, concluding re¬ 
marks are given in Section^ 


2. Prominence thread model and numerical solution 
of the eigenvalue problem 

The equilibrium model of a prominence thread is schematically 
shown in Figure Q] We use a cylindrical coordinate system, with 
r, ip, and z the radial, azimuthal, and longitudinal coordinates, re¬ 
spectively. We consider a cylindrically symmetric straight mag¬ 
netic flux tube of radius R and length L. The ends of the tube are 
located at z = ±L/2 and are line-tied at two rigid walls represent¬ 
ing the solar photosphere. The centre of the tube corresponds to 
z = 0. The magnetic field, B, is straight and along the axis of the 
tube. The magnetic field strength, /i (l , is constant everywhere. 
We use the f3 - 0 approximation, where f3 refers to the ratio of 
the gas pressure to the magnetic pressure. This is an appropri¬ 
ate approximation to study kink oscillations^ which are mainly 
driven by magnetic tension in thin tubes (iGoossens et al.ll2009i) . 
In the p — 0 approximation we can freely choose the spatial dis¬ 
tribution of mass density. Hence, the equilibrium density, po, is 


Po(r,z) 


Pi(z), if r < R, 
Pe(z), if r > R, 


( 1 ) 


where the internal, pi(z), and external, p e (z), densities are func¬ 
tions of z only, and p,fz) > p e (z). The external plasma rep¬ 
resents the coronal medium. In the present model the density 
jumps abruptly at the boundary of the tube, i.e., at r = R. We 
are aware that a continuous transition of density in the radial di¬ 
rection would be more realistic, although it would complicate 
matters. The kink oscillations of a transversely nonuniform tube 
are resonantly coupled to Alfven waves and, consequently, the 
oscillations are damped in time (see, e.g.. Goossens et al . ?snb. 
In such a case, the kink mode is no longer a true normal mode 
of the flux tube but a spectral pole of the Green’s function, re¬ 
lated to the initial-value problem, that represents a damped co¬ 
ordinated motion of the plasma ( ISedlaceklll97ll ). Since in the 
present work we are not interested in the damping, we assume a 
jump of density at the boundary of the tube to avoid the resonant 
absor ption process. Readers are ref erred to lArregui et al.l 12008. 
1201 ll) and ISoler et ail (120091 1201 Oh for more information about 
the resonant damping of prominence thread kink oscillations and 
to IGoossens et al.l (1201 ll) for a general review. A comparison 
of the efficiency of several damping_mechanisms of transverse 
thread oscillations is given in ISoler et al.l (|2014|) . 

We denote the internal density at the centre and at the end 
of the flux tube as pi(0) = pj,o and p;(L/2) = pix/ 2 , respectively. 
Then, we define the ratio of the two densities as* = Pi.o/Pi.L/ 2 , 
with * > 1. Contrary to coronal loops, prominence threads are 
denser near the centre of the thread than at their footpoints. The 
larger*, the stronger the density variation along the thread, with 
the case * = 1 corresponding to a homogeneous thread. Here, 
we consider three different longitudinal profiles for the internal 
density, namely a Lorentzian profile 


, \ _ _ Pio _ 

P ' ’ l+4(x-l )z 2 /F 2 ’ 

a Gaussian profile 


Pi(z) = pi, 0 exp 


-4 (log*) 



and a parabolic profile 

Pi(z) = Pi,o 11 


( 2 ) 


(3) 


(4) 
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The three density profiles are compared in Figure [2] The three 
profiles represent different plasma arrangements within the mag¬ 
netic tube. In the Lorentzian profile the densest plasma is nar¬ 
rowly concentrated near the centre of the tube, while the rest of 
the tube is occupied by much less dense plasma. In essence, the 
Lorentzian profile i s similar to the p i ecewi s e models used in pre- 
vious works (e.g., iDiaz et aU 120021120101: iTerradas et all [2008: 
ISoler et aDl2010l) . On the contrary, in the parabolic profile the 
dense prominence plasma is broadly distributed along the tube. 
The Gaussian profile represent an intermediate situation between 
the Lorentzian and parabolic profiles. The actual density varia¬ 
tion along prominence threads is probably linked to the process 
that leads to the formation of prominences and the condensation 
of the cool and dense plasma. A thermal instability might be in¬ 
volved in that process (see, e.g., [Luna et akll2012l) . The purpose 
for choosing these three paradigmatic profiles is to determine 
the effect of different plasma arrangements on the period ratio of 
kink oscillations. 



Fig. 2. Variation of the internal density along the prominence thread for 
the Lorentzian profile (black solid line), the Gaussian profile (red dotted 
line), and the parabolic profile (blue dashed line). We used^ = 100 for 
representation purposes. 


On the other hand, since we focus on horizontal threads, 
the external coronal plasma is assumed uniform for simplicity, 
namely p e (z) = p e . Then, we define the prominence to corona 
density contrast as the ratio of the internal density at the centre 
of the thread to the external density, namely £ = pio/p e . 


The linear ideal MHD equations for a /? = 0 plasma that 
govern the behavior of small perturbations superimposed on the 
equilibrium prominence thread are 

= -(Vxb)xB, (5) 

F 

= V x (v x B), (6) 

where v = (v r , v v , v z ) and b = (b r , b z ) are the velocity and 
magnetic field perturbations, respectively, and p is the magnetic 
permeability. Due to the f3 = 0 approximation, v, = 0. Since the 
equilibrium is invariant in the azimuthal direction, we can put the 
perturbations proportional to exp ( irrnp ), where m is the azimuthal 
wavenumber. We consider kink oscillations, so we set m = 1. In 
addition, we perform a normal mode analysis and express the 
temporal dependence of the perturbations as exp (-icot), where 
u> is the frequency. The period of the oscillation is P = 2njo>. 
Then, Equations (0 and 0 become 


d\ 

P °di 

db 

dt 


(jJV r 

m )' 

( dz 

* 

OJVy 

Bo | 

m>' 

ft; 

\r " 

Lob* r 

_ R 9v r 

tJby 

_ -ft* 

- B0 'AT’ 

cob: 

/ dv r 


d_K\ 

dr )' 

dbA 
dz /’ 



(7) 

( 8 ) 
(9) 

( 10 ) 

(ID 


where, with no loss of generality, we defined the new perturba¬ 
tions v* = ivy, b* = ib r , and b* = ib z in order to remove complex 
numbers and so to have equations with real coefficients. This is 
computationally convenient. Physically, the factor i accounts for 
a phase difference of n/2. Therefore, the new perturbations v*, 
b*, and h*. are 7r/2 out of phase with respect to the actual pertur¬ 
bations v v , b r , and b z . The remaining perturbations, namely v r 
and b v , are not altered. From here on, we drop the superscript * 
for simplicity. Equations (I71i (fTH define a 2D generalized eigen¬ 
value problem, where co is the eigenvalue and the perturbations 
form the eigenvector. Contrary to the piecewise constant mod- 
els fr e quently used in the literature (e.g.jDvmov a & kud ermanl 
120051: ISoler et al.l|2010t iDfaz et all 120101) . the analytic solution 
to the eigenvalue problem becomes very c omp licated when po 
is an arbitrary function of r and z (see lAndries et al..'2005h ). In 
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this work, the 2D eigenvalue problem is numerically solved with 
the PDE2D code dSewellll2005l) . a general-purpose partial differ¬ 
ential equation solver. The PDE2D code was_previously used to 
solve a similar 2D eigenvalue problem by [Arrestui et al . ( 201D . 
The PDE2D code uses a collocation method, and the general¬ 
ized matrix eigenvalue problem is solved using the shifted in¬ 
verse power method. The code uses finite elements and allows 
the use of non-uniformly distributed grids. Different grid resolu¬ 
tions have been tested so as to assure the proper convergence of 
the solutions. The code provides the closest eigenvalue to an ini¬ 
tially provided guess and the corresponding spatial form of the 
eigenfunctions. The boundary conditions used in the code are 
those consistent with trapped and standing kink oscillations. The 
two non-zero components of the velocity perturbation, namely 
v r and v v , and the longitudinal component of the magnetic field 
perturbation, namely b z , vanish at z = ±L /2 because of the line- 
tying boundary condition at the photosphere. In turn, b r and b^ 
have vanishing longitudinal derivatives at z = ±L/ 2. In the radial 
direction, v r , v v , b r and have vanishing radial derivatives and 
b z = 0 at the axis of the tube, r = 0. The code self-consistently 
connects the internal perturbations with the external perturba¬ 
tions at r = R and satisfies the jump relations that naturally arise 
from the integration of the equations across r = R. The con¬ 
dition that the modes are trapped imposes that all perturbations 
vanish far away from the tube in the radial direction, i.e., when 
r —> oo. This condition is accomplished in the numerical code 
by setting all perturbations equal to zero at r — r max , where r max 
is the largest value of the radial coordinate considered in the nu¬ 
merical domain. Ideally, we should set r max —> oo. However, this 
is not possible computationally, and r max must be set to a finite 
value. We locate r max sufficiently far from the tube boundary to 
properly compute the drop-off rate of perturbations in the radial 
direction and so avoid numerical errors. To make sure that r max 
is sufficiently large, we have performed convergence tests by in¬ 
creasing r max until the solutions showed no dependence on this 
parameter. We took r max = 20 R in the computations we include 
in this paper. In summary, the boundary conditions used in the 
numerical code are 


= 0, b z = 0, 


dv r dv<p db r db^ 

d7 ~ ~dr ~ ~d7 “ Hr 

v r = = by = b tp = b z = 0, 

db r db^ 


v r = Vu> = b z = 0, 


dz dz 


= 0 , 


at 

r = 0, 

(12) 

at 

^ — y max? 

(13) 

at 

z = +L/2. 

(14) 


In the specific case that both the internal and external densi¬ 
ties are uniform, namely pi(z ) = p\ and p e (z) = p e , it is possi¬ 
ble to derive analytic expressions for the period of kink oscilla- 
tions in the thin tube (TT) limit, i.e., when L/R » 1 (see, e.g., 
[Edwin & Robertsll983fclGoossens et al,l2009h . The period of the 
fundamental longitudinal mode, Pq , and that of the first longitu¬ 
dinal overtone. Pi, are given by 


. = 2l un p = a. im 

v A ,i y 2£ ’ v A ,i yj 2£ 


(15) 


where v A ,i = Bo/ yfjlpi is the internal Alfven velocity. Thus, 
the period ratio is Pq/P\ = 2 in longitudinally homogeneous 
thin tubes. In the present model, deviations from Po/P\ = 2 
can be caused by dispersion when the thickness of the thread 
does not satisfy the TT condition and by longitudinal varia¬ 
tion of density. The departure from the TT limit is not likely 
to have a strong impact for prominence threads, since observa¬ 
tions clearly show that threads are very thin and long structures 


(see lLinlf201 1). Although only the part of the tubes filled with 
the densest and coolest plasma can be seen in, e.g.. Her and Ca II 
observations, the length of the whole magnetic tube, L, must be 
much longer than the observed length of the prominence threads, 
Lnhs- The values of R and L„ h s reported by the o bservations (e.g., 
Oka moto et al .120071 iLin et all2008tlLinil201 lh are in the ranges 
50 km < R < 3 00 km and 3000 km < L 0 b s < 28,000 km, re¬ 
spectively, while lOkamoto et all d2007l) estimated L to be at least 
L ~ 10 5 km. These numbers give 10 < L a o s /R < 400 and 
400 < L/R < 2000, meaning that real prominence threads are 
actually very thin tubes. Therefore, in the remainder of this pa¬ 
per we shall focus on the effect of longitudinal density variation, 
which may have a much stronger impact on the period ratio. 


3. Results 

We start by investigating the modification of the fundamental 
mode and first overtone periods when the ratio of the central 
density to footpoint density, x, is increased from x = I (ho¬ 
mogeneous tube) to^- = 100. This is achieved by reducing the 
density at the footpoints while the central density is kept fixed. 
These results are displayed in Figure[3] We denote by P' and P' 
the fundamental mode and first overtone periods, respectively, 
numerically computed with the PDE2D code. In turn, Pq and P\ 
are the periods of a homogeneous tube with p\ = p\ o computed 
from Equation (031 ). 

When x increases, both the fundamental mode (Figure [3}i) 
and first overtone (Figure [3 Jd) periods decrease with respect to 
the values for a homogeneous tube with a density equal to the 
central density since, in average, the internal density is reduced. 
The decrease of the periods depends on the specific density pro¬ 
file considered. In this regard, the parabolic profile produces the 
smallest decrease of the periods. The results for the parabolic 
profile saturate at x ~ 10 and become independent of x af- 
terwards. On the contrary, the Lorentzian profile produces the 
largest decrease of the periods, which keep decreasing as x in¬ 
creases. The results for the Gaussian profile are in between those 
of the Lorentzian and parabolic profiles. 

Figure [3J: shows the dependence with x of the numerically 
computed period ratio P/J P \. Consistently, P' } / P' = 2 when 
X - 1. As ^ increases, P'q/P\ takes values larger than 2. This re¬ 
sult conceptually agrees with that obtained bv lDIaz et alJ ( 2010 ) 
in tubes with a piecewise constant density profile. Here, we find 
that Pg/Pj > 2 in the case of continuous density profiles as well. 
Again, how large is the increase of the period ratio from 2 de¬ 
pends strongly on the density profile considered. The Lorentzian 
profile produces a significant increase of the period ratio, so that 
Pg/Pj » 3.5 when x = 100. Conversely, the increase of the 
period ratio is small for the parabolic profile and P' /P' ~ 2.2 
when;y = 100. Finally, we see again that the results for the Gaus¬ 
sian profile are in between those of the Lorentzian and parabolic 
profiles, so that the period ratio increases moderately with x for 
Gaussian profile and is P' ( /Pj ~ 2.5 when^ = 100. 

It should be noted that comparing the influence of the various 
profiles using the results of Figure[2is not straightforward. When 
the ratio of the central density to footpoint density is fixed, the 
total mass in the thread is different for every profile. This fact 
makes a direct comparison between profiles unsuitable. To fairly 
compare the results for the various profiles, the total mass in the 
thread should be the same, which is equivalent to considering 
the same average density in the thread (see lAndries et al.l2005al) . 
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Fig. 3. Dependence of the numerically computed (a) fundamental mode 
period, P' 0 , (b) first overtone period, P \, and (c) period ratio, , P' 0 /P\, 
with x f° r the Lorentzian profile (black solid line), the Gaussian profile 
(red dotted line), and the parabolic profile (blue dashed line). We used 
L/R = 100 and f = 100 in all cases. P 0 and Pi denote the periods of a 
homogeneous thread with p\ = pi.o. 


Fig. 4. Same results as Figure [3] but as function of the ratio of the av¬ 
erage internal density to the central density, ( p\) /pyj. For comparison, 
the green dashed-dotted line shows the result for a homogeneous thread 
with pi = (pi), i.e., the ‘average thread’. The symbols in panel (c) cor¬ 
respond to the empirical fit of Equation ( 1171 . 


Fience, we define the average internal density as 

i r L/2 

(Pi) = — Pi(z)dz. (16) 

L J-L/2 

Figure |4]displays the same results as Figure [3] but now (p{) /pm 
replaces x hi the horizontal axes of the plots. We see that there 


are no significant differences between the curves for the various 
profiles when the average internal density is the same. Impor¬ 
tantly, Figure[4j: shows that the behavior of the period ratio is the 
same for the three profiles. The period ratio increases as (p,) /p,o 
decreases. This result indicates that the ratio of the average in¬ 
ternal density to the central density, namely (pj) /p K o, and not the 
ratio of the central density to the footpoint density, namely x > is 
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Fig. 5. Average internal density, {p\) /pio, as a function of^. The vari¬ 
ous lines are for the Lorentzian profile (black solid line), the Gaussian 
profile (red dotted line), and the parabolic profile (blue dashed line). 


the parameter that actually controls the behavior of the period 
ratio regardless of the form of the density spatial variation along 
the thread (see lAndries et al.ll2005ai . for an equivalent result in 
the case of coronal loops). We performed an empirical fit to the 
results of Figure [4]e and obtained that the dependence of the pe¬ 
riod ratio on the average density is very closely approximated by 


li + (Mf 2 

p\ Uo/ 


(17) 


for all the profiles. The empirical Equation (IT71) is overploted in 
Figure[4j: to show that the fit is excellent. Equation (IT71) is inde¬ 
pendent of the density profile. It is an important result that has 
simple and straightforward applications for prominence seismol¬ 
ogy. The implications of this result for prominence seismology 
are discussed in Section[4] 

Note that the length of the curves in Figure |4] is different 
for different density profiles. The reason for this is that (p;) was 
computed by varying^- in the interval [1,100], but the depen¬ 
dence of the average density on x is different for every profile. 
The analytic result of computing the average density from Equa¬ 
tion © is 


arctan yj x - 1 
<Pi> = Pi,o- 

for the Lorentzian profile. 



<A> = Pi,0 


Vfr erf ^ /logx 
2 V lo 8T 


(18) 


(19) 


for the Gaussian profile (where erf and log denote the error func¬ 
tion and the natural logarithm, respectively), and 


<Pi> = Pi.O 


2y+l 

3t 


( 20 ) 


for the parabolic profile. These average densities are plotted in 
Figure[5]as function of We see that the average density of the 
Lorentzian profile is lower than that of the parabolic profile for 
the same value ofyo with that of the Gaussian profile in between. 
For prominence threads it is expected x to be a large parameter. 


Hence, we perform the limit x » 1 and the expressions for the 
average density of the different profiles simplify to 


<Pi> ~Pi,o|* 1/2 , (21) 

for the Lorentzian profile, 

<Pi> ~ Pi.o-^ (logT) 1/2 . (22) 

for the Gaussian profile, and 

<Pi>~Pi,o^, (23) 

for the parabolic profile. Importantly, the average density of the 
parabolic profile becomes independent of x when x » 1 ■ Now 
we substitute these average densities into Equation (ITTI ) to obtain 
simple relations between P'JP\ and x for the various profiles, 
namely 



for the Lorentzian profile. 


n 

p\ 


~ i + 



for the Gaussian profile, and 


(24) 


(25) 



(26) 


for the parabolic profile. These approximate results agree well 
with the numerical solutions displayed in Figure [3] when x » 1. 

For comparison, we overplot in Figure[4]the results for a lon¬ 
gitudinally homogeneous thread whose internal density is equal 
to the average density, i.e., p\ = (p;). We call this hypotheti¬ 
cal case the ‘average thread’. Obviously, the period ratio of the 
‘average thread’ is 2 due to the absence of longitudinal inho¬ 
mogeneity. However, we see that the fundamental mode period 
of the ‘average thread’ is shorter than the period of a nonuni- 
form thread, but the first overtone periods are similar in both 
cases. The fact that the density varies spatially is important for 
the fundamental mode period, but is not very relevant for the 
first overtone period. This points out that ignoring the longitu¬ 
dinal variation of density in prominence threads by using ‘av¬ 
erage thread’ models may provide inaccurate values of the fun¬ 
damental mode perio d. This result is op posite to that found by 
lAndries et al] d2005al) for coronal loops. lAndries et al] (|2005a) 
showed that the first overtone was more affected by plasma strat¬ 
ification than the fundamental mode. The reason for this differ¬ 
ence is probably the fact that in prominence threads the densest 
plasma is located near the tube centre, while in stratified coronal 
loops the pl asma is denser near the footpoints. Also note that, 
contrary to lAndries et ab ( 2005a ). we do not perform a weighted 
spatial average of the density but an unweighed average. 

Finally, for completeness we have computed the spatial form 
of the eigenfunctions with the PDE2D code. We display in Fig¬ 
ure [6] longitudinal cuts of v r , and b z near the boundary of 
the thread, i.e., r * R, for both the fundamental mode and the 
first overtone using a particular set of parameters. Figure[6]may 
be compared with the eigenfunctions in piecewise-like models 
(ISoler et alJl2010t lArregui et al]|201 lh . We can observe in Fig¬ 
ure [6] that the eigenfunctions of the fundamental longitudinal 
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Fig. 6. Longitudinal cuts at r = R of v r (top panels), v v (mid panels), and b- (bottom panels) corresponding to the fundamental mode (left panels) 
and the first overtone (right panels). The various lines are for the Lorentzian profile (black solid line), the Gaussian profile (red dotted line), and 
the parabolic profile (blue dashed line). We used L/R = 100, x = 50, and = 100 in all cases. Arbitrary units are used. 


mode of transverse kink waves (left panels) have no nodes in 
the axial (longitudinal) direction, while the eigenfunctions of 
the first longitudinal overtone of transverse kink waves (right 
panels) have one node at the centre of the tube. We find that 
v r and v<p corresponding to the different profiles are very simi¬ 
lar for both fundamental mode and first overtone. This suggests 
that measurements of the thread transverse velocity or, equiv¬ 
alent, the displacement could provide little information about 
the density profile. However, the spatial form of b z is sensitive 


to the density profile. The more concentrated the dense promi¬ 
nence plasma near the centre of the tube, the faster the ampli¬ 
tude of b : drops from the centre. This result was also noticed by 
lArregui et al.l (1201 ll see their Figure 6c). In the present compu¬ 
tations, b z is essentially confined near the centre of the thread for 
the Lorentzian profile, but it is almost a sinusoidal-type function 
for the parabolic profile. 
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4. Prominence seismology 

The results of this paper and, in particular, the empirical fit 
of Equation (IT71) have important implications for seismology 
of prominences. Equation (TT71 ) is general and does not depend 
on the specific density profile. This means that it can be used 
even when the density profile is not known. By inverting Equa¬ 
tion o. we find that the ratio of the average internal density to 
the central density is related to the period ratio as 


P' 

^-1 
P, 


-2 


(27) 


(Pi) 

Pi,0 \- i 

Equation (l27l > can be used to seismologically estimate the ratio 
(pi) /pi o in prominence threads using observations of the period 
ratio. 

As mentioned in the Introduction, no reliable simultaneous 
observations of the two periods in prominence threads are cur¬ 
rently available. Two periods were reported bv lLin et~ai] (l2007i) . 
These authors studied an extensive area in a prominence and it is 
unclear whether the two periods originate from the same thread 
or correspond to different threads. However, to perform a simple 
exemplifying application, let us proceed as i nlDIaz et ah] (2010b 
and assume that the two periods reported bv lLin et al.l ( 2007 ) ac¬ 
tually correspond to the fundamental mode and the first overtone 
of the same oscillating thread. The two periods are P' 0 = 16 min 
and P\ = 3.6 min, which give P' 0 /P[ ~ 4.44. Then, using Equa¬ 
tion (1271) we get (p;) /p; o ~ 0.084. This result indicates that there 
is a strong density gradient along the thread, since the average 
density is a small fraction of the central density. This is compat¬ 
ible with a large value of the ratio of the central density to the 
footpoint density, i.e.,^ » 1. 

Next, we can also use the observed P' Q /P\ to infer the ratio 
of the central density to the footpoint density, x- To do so, we 
must necessarily assume a particular density profile along the 
thread. Note that the form of the density profile was not needed 
to estimate (p;) /p; o- In the case of the Lorentzian profile we get 
from Equation (l24l) that 


X : 




4 \P\ 


347, 


(28) 


while for the Gaussian profile we get from Equation (1251) that 


X ~ exp 


-i5-i 


4 \P' 


10 


48 


(29) 


The value of x corresponding to the Gaussian profile is unreal¬ 
istically large, while the value corresponding to the Lorentzian 
profile is reasonably consistent with the expected densities in 
prominence threads. In the case of the parabolic profile, we 
found in Equation (l26l > that its period ratio is 2.22 when x » 1, 
so that the observed period ratio of 4.44 seems hardly compati¬ 
ble with a parabolic profile. The inferred values ofx suggest that, 
among the three profiles studied here, the Lorentzian profile may 
provide the better explanation for the ratio of the two periods 
reported by iLin et all ( 2007 ). However, a rigorous comparison 
between models should be done and an adequate treatment of 
errors should be considered to obtain more robust conclusions 
(see lArregui et al1l2013l) . 

5. Concluding remarks 

In this paper, we have investigated the ratio of the period of the 
fundamental mode to period of the first overtone of kink os¬ 
cillations of longitudinally inhomogeneous prominence threads. 


Contrary to previous works that used piecewise constant den¬ 
sity profiles (iDIaz et al.ll2010i) . here we assumed a continuous 
nonuniform density profile along the thread, that is more consis¬ 
tent with the simulations of thread formation viajplasma conden¬ 
sation in solar prominences (e.g., lLuna et al.ll2012b . Our results 
indicate that the ratio of the average internal density to the cen¬ 
tral density, namely (pi) /pi.o, is the parameter that controls the 
behavior of the period ratio regardless of the form of the density 
spatial variation along the thread. A fit of the period ratio as a 
function of (pj) /p K o was obtained, which can be used to estimate 
the degree of nonuniformity along the threads if observations of 
the two periods are reported. From the observational point of 
view, the detection of the first overtone is challenging, since the 
amplitude in the densest part of the tube is small for that mode. 
Although difficult, we believe that the observation could be done 
with present or future high-resolution instruments. For instance, 
using recent observations with the Hinode satellite, iHillier et ali 
(12013l) reported a large number of oscillations in vertical promi¬ 
nence threads with velocity amplitudes as low as 0.2 km s' 1 . 
In view of these observations, the detection of oscillations with 
such small amplitudes in horizontal threads is more than likely. 

Here we have studied the period ratio from a purely nu¬ 
merical point of view. Hence, the equation that relates the pe¬ 
riod ratio with the average density (Equation (ITTI O was empir¬ 
ically obtained from the numerical data. Although it is beyond 
the purpose of the present study, it would be desirable that fu¬ 
ture works pursue a more robust theoretical, i.e., analytical, ex¬ 
planation for this relation. The analytical theory developed by 
iDvmova & Rudermanl (120051) could be used to do so. 

In the model, we assumed for simplicity that the density 
within the flux tube is symmetric with respect to z = 0. The ef¬ 
fect of shifting the density enhancement from the tube centre was 
investigated in a previous paper bv ISoleretak (2 01 Oh using a 
simple piecewise model. The conclusion from lSoler et al.l ( 2010 ) 
was that the periods of kink oscillations are weakly affected by 
the position of the densest plasma within the tube unless the den¬ 
sity enhancement is put very close to the footpoints, which is not 
realistic in the case of horizontal prominence threads. Therefore, 
we expect that the use of an asymmetric density profile would 
not change significantly the results of the present paper. 

The results of the seismology application of Section [4] sug¬ 
gest that value of the period ratio could be used to distinguish be¬ 
tween different density profiles. In a general sense, rigorous sta¬ 
tistical methods to perform model comparison, as the Bayesian 
analysis used by [Arregui et al.l d2013l) . could be very useful to 
determine what kind of density variations are the most probable 
ones for a given value of the period ratio. Future works should 
exploit this possibility. Also, the presence of uncertainties in the 
measurements of the periods should be properly accounted for. 

Finally, we should mention that there are some effects 
not included in the prominence thread model that may in¬ 
fluence the value of the period ratio. For instance, an in¬ 
gredient missing from the model is the presence of mass 
flows along the thread. Simultaneous transverse oscillations 
and flows in prominence fine structures have been observed 
(e.g., iOkamoto et al.l 12007) . The effect of flows on the pe¬ 
riod ratio has been investigated in piecewise cons tant mod¬ 
els (e.g., ISoler & Goossensll201 lfc lErdelvi et al.ll2014h . but not 
in threads with nonuniform density profiles. Among other ef¬ 
fects, the roles of magnetic field twist, curvature, and cross sec- 
tion have been studied in the case of coronal loops (see, e.g. , 
Verth_&Erdel^i 2008; Vertheta^ 2008|: Ruderma n etal.l feoOS: 


Morton & Erdel vil2009( Karami & Baharil2012l) but they should 

be explored for prominence threads as well. In addition, in view 
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of the large amount of observations (see, e.g., iHillier et alJfoTl . 
the problem of the transverse oscillations of vertical threads is 
also theoretically relevant. For vertical threads, the effect of grav¬ 
itational stratification should be included in the model. A com¬ 
parison between the properties of the oscillations of vertical and 
horizontal threads is an interesting task to carry out in forthcom¬ 
ing works. 
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Appendix A: Comment on the paper by 
Lomineishvili et al. (2014) 

The results computed in the present paper for a thread with 
a paraboli c density profile are in disa greement with those ob- 
tained by ILomineishvili e t al.l (1201 4l) using the same model. 
Lomineis hvili et alJ (12014 ) claim that the period ratio is 3, but 
we found that the correct value of the period ratio is 2.22, ap¬ 
proximately. Here we explain the source of the discrepancy be¬ 
tween the results of lLomineishvili etal . d2014l) and those of the 
present paper. 

A property of the analysis by ILomineishvili et all (12014l) is 
the discontinuous behaviour of their solutions as a function of 
the inhomogeneity parameter a at a — 0. The parameter a used 
by ILomineishvili et akl (12014l) is related to our parameter x as 
a - (x - 1) !x, so that a = 0 corre sponds to x - 1 in our nota¬ 
tion. As ILomineishvili et al J (20141) observe (see the paragraph 
that follows their Equation (19)) a = 0 in their Equation (17) 
leads to the solution of the homogeneous tube. They continue 
to say that after corresponding calculation, one may recover the 
well-k nown dispersion relation obtained by lEdwin & Robertsl 
( 198 31). The ratio of the period of the fundamental longitudinal 
mode to that of the first overtone for a homogeneous tube in the 
T T approximation is 2 . For a + 0, however small the value of a 
is. ILomineishvili et alJ d2014l) claim that thi s ratio is 3 in the TT 
approx imation. According to the analysis bv lLomineishvili et alJ 
d2014l) the period ratio varies in a discontinuous manner at a — 0, 
where it jumps from 2 for a - 0 to 3 for a > 0. Since longitu¬ 
dinal density stratification is not a singular perturbation to the 
problem of linear waves on magnetic cylinders, this result is in¬ 
consistent. 

Equations (53) and (54) of lLomineishvili et al.l ( 2014 ) corre¬ 
spond to analytic expressions for the fundamental mode and first 
overtone frequencies in the TT approximation. Setting a —» 0 in 
Equations (53) and (54) of lLomineishvili et al.l ( 2014 ) results in 
vanishing frequencies, i.e., infinite periods. This result can also 
be seen i n the numerically obtained frequencies plotted in their 
Figure 3.ILomineishvili et al. j|2014l) do not recover the frequen¬ 
cies of lEdwin & Robertsl d 19831) when a —> 0. 

According to | Z aqarashviIil d2014t) . the modes obtained by 
ILomineishvili et all (l2014i) 'owe their existence to inhomogene¬ 
ity and so they disappear when the thread is homogeneous, 
i.e., when a —> 0. Also according to Zaqa rashvili d2014l) . 
ILomineishvili et all d2014l) assumed open boundary conditions 
at the ends of the thread because in their model “the thread it¬ 
self probably is a part of much longer magnetic tube, which has 


a coronal density outside the thread”. However, fixed boundary 
conditions are necessary fo r standing kink oscillation s. Hence, 
the mathematical analysis of lLomineishvili et al.l d2014l) does not 
lead to a dispersion relation valid for standing kink oscillations. 
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